Quantum critical temperature of a modulated oscillator 
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We show that the rate of switching between the vibrational states of a modulated nonlinear 
oscillator is characterized by a quantum critical temperature Tc oc h^ . The rate is independent 
of T for T < Tc- Above T^ there emerges a quantum crossover region where the slope of the 
logarithm of the distribution over the oscillator states displays a kink and the switching rate has the 
Arrhenius form with the activation energy independent of the modulation. The results demonstrate 
the limitations of the real-time instanton theory of switching in systems lacking detailed balance. 



PACS numbers: 05.40.-a, 03.65.Yz, 74.50.+r, 85.25.Cp 

I. INTRODUCTION 

Thermal equilibrium systems have time reversal sym- 
metry, which leads to detailed balance: the rates of tran- 
sitions back and forth between any two states are equal 
[T]. Nonequilibrium systems generally do not have de- 
tailed balance. An important exception is a nonlinear os- 
cillator additively modulated close to its eigenfrequency 
LjQ or parametrically modulated at ~ 2ujo [3]. This 
system has detailed balance in the rotating wave approx- 
imation (RWA), but only for zero temperature [1]. For 
nonzero temperature the oscillator does not have detailed 
balance, cf. [5]. 

Interestingly, the detailed-balance probability distribu- 
tion of the oscillator is "fragile". The distribution for 
T — > can be exponentially strongly different from that 
for T = for weak damping and comparatively strong 
modulation, where the oscillator dynamics is semiclassi- 
cal [51 [7] . The fragility was found in the region where the 
oscillator has two stable vibrational states (SVSs). As 
its consequence, the rate of switching from a vibrational 
state Wsw can be exponentially increased for T > com- 
pared to its value for T = 0. Such increase is of significant 
interest, as the bistability of forced vibrations of quan- 
tum nonlinear oscillators plays an important role in their 
applications in quantum information, and the interstate 
switching is attracting much attention, cf. [5Hl3). 

In this paper we show that the T = region and the 
region, where the lack of detailed balance in the under- 
damped modulated oscillator is pronounced, are sepa- 
rated by a quantum temperature oc fP. This temper- 
ature does not show up in the standard WKB approxi- 
mation, it is not related to the conventional WKB cor- 
rections oc h^. For T > Tc the exponent of the switching 
rate Wgw changes with increasing T from the semiclassi- 
cal T = to the semiclassical T value. The change 
occurs in the temperature range that scales with h as 

h/\\ogh\. 

The distribution of the modulated oscillator over its 
quantum states is formed as a result of the coupling to a 



thermal reservoir. The coupling leads to oscillator relax- 
ation. In a simple picture relaxation comes from emis- 
sion of excitations of the bath, for example, photons. It 
is invariably accompanied by noise, because photons are 
emitted at random. For T > 0, along with the noise 
from photon emission, there is noise from photon ab- 
sorption. Its relative intensity for low T is oc n, where 
n = [exp{hu!o/kBT) — is the oscillator Planck num- 
ber. The noise leads to fluctuations of the oscillator and 
ultimately to switching between the SVSs over an effec- 
tive barrier in phase space [51 [7] . 




FIG. 1. (a) The scaled Hamiltonian g{Q,P) of the reso- 
nantly driven nonlinear oscillator; Q and P are the dimen- 
sionless coordinate and momentum in the rotating frame. In 
the presence of weak dissipation, the local maximum and min- 
imum of g correspond to the small- and large-amplitude sta- 
ble vibrational states (SVSs). The scaled field intensity is 
/3 = 0.02. (b) The cross-section g{Q,0); the solid lines show 
the quasienergy states localized about the local maximum of g 
at Q = Qsvs f'^'' ^ ~ 0.015. The scaled quasienergies of these 
states gn lie between pg^g = g{Qsvs'^) ^^'^ 9s ~ ^(QsiO), 
where {Q^ ,P = Q) is the saddle point of g{Q, P). 

An insight into the onset of the switching rate fragility 
can be gained from Fig. [l] which shows the effective os- 
cillator Hamiltonian in the rotating frame. Its extrema 
correspond to the SVSs, for weak damping. Quantum 
states of the oscillator localized about an extremum of 
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the Hamiltonian are analogous to intrawell states of a 
particle in a static potential. For weak damping, the 
state populations pn are determined by a balance equa- 
tion. In dimensionless time r, see Eq. ([5]), it reads 

Pn = dpn/dr = '^{WmnPra " WnraPn)- (1) 
m 

The dimensionless transitions rates Wmn = + 
Wmif^ have two terms, which come from emission and 
absorption of bath excitations, respectively; Wmrf^ oc n. 

For T = 0, the stationary probability distribution pi"'* 
is formed by the rates Wmn- Of interest for switching is 
the population of quantum states n with the wave func- 
tions extending to the boundary of the classical basin 
of attraction to the SVS. For such states n ^ 1, for small 
h (we count the states off from the SVS). In a certain 
parameter range the rates Il^m n'''' of absorption-induced 
transitions decay with n — m > slower than pn^ / p^m ■ 
Then, if we consider these transitions as a perturba- 
tion in Eq. the correction Y^rn<n^rnn^^ Pm / Pn'' to 

Pn/ Pn diverges for n — >■ oo. This correction describes 
the absorption-induced flux into state n from the states 
that are closer to the SVS. The divergence occurs even 
for n — >■ 0, indicating the fragility of the T = solution. 

The value of Tc can be estimated by noticing that, for 
finite fi., the number M of quantum states localized about 
an SVS is finite, M oc Therefore the above correc- 
tion has a finite number of terms. If the characteristic 
ratio of the consecutive terms is w > 1, the correction is 
c>c nw*^ for n ~ M. It remains small for small nw*^, i.e., 
for T <Tc^ {hujo / k b)/M log V cx ff. 

Finding the distribution for a finite number of states 
requires calculating the rates Wnm of transitions between 
remote states, \n — m\ ^ 1. We do this for an oscillator 
modulated close to its eigenfrequency by combining the 
WKB and the conformal mapping techniques. 

II. OSCILLATOR DYNAMICS IN THE 
ROTATING FRAME 

The Hamiltonian of the modulated oscillator reads 

Ha = Y + l^oQ^ + ^q^-Aq cosiupt), (2) 

where q and p are the oscillator coordinate and mo- 
mentum. We assume that the modulation detuning 
doj = u)F — uq is small, \Su;\ <^ loq, as is also the non- 
linearity , |7|((i'^) ^ The oscillator displays bistabil- 
ity of forced vibrations for ^ 5uj > 0; for concreteness we 
assume 7, Slu > 0. 

We switch to the rotating frame using the transforma- 
tion U{t) = exp{—iujpa^ at), where and a are the oscil- 
lator raising and lowering operators. We also introduce 
slowly varying in time dimensionless coordinate Q and 
momentum P, W[t)qU{t) = C (Q 



and W{t)pU{t) = —Cujf {Q sinujpt — P cosutpt) , with 
the scaling constant C — {Su}p6uj/3"fy^^; 

[Q, P] = zA, A = 3h-f/8Lul{ujp - Luo). (3) 

The dimensionless parameter A plays the role of the 
Planck constant in the dynamics in the rotating frame. 
We assume A ^ 1, so that this dynamics is semiclassical. 

In the RWA the transformed oscillator Hamiltonian 
Ho = [8w|.((5aj)^/37].g is independent of time. Here, 

g^9{Q,P)^\{Q^ + P'-lf--f3'/'Q; (4) 

13 = SIjIA^ /32ujp\ujp — Wop is the scaled field intensity. 
Function g{Q,P) is shown in Fig.[l] Where the oscillator 
is bistable, g{Q, P) has the form of a tilted Mexican hat 
with the extrema corresponding to the SVSs. 

Operator g is Hermitian and has a complete set of 
eigenfunctions ipniQ), gi^n = gn'4'n- The eigenvalues (?„ 
give the scaled quasienergies £„ of the modulated oscilla- 
tor in the RWA, e„ sa [8ci;|.((5a;)^/37](7„. They are shown 
in Fig. [l|b) . The spacing between the eigenvalues (?„ is 
cx A. The number of states between an extremum and 
the saddle point of g{Q, P) is M ^ A~^ > 1. 

A. Master equation and the transition rates 

The dynamics of the oscillator weakly coupled to a 
thermal reservoir can often be described in slow time by a 
Markov master equation for the oscillator density matrix 
p. To the lowest order in the coupling 

dp/dr = -iX^^ [g, p] ~ kp, t = \ujp - uja\t. (5) 

Operator kp describes relaxation; the coupling-induced 
renormalization of the oscillator parameters is incorpo- 
rated into Wo, 7. If the coupling is linear in the oscillator 
raising and lowering operators, kp — K{n + l){a^ap — 
2apa^ + pa^a) + Kn{aa) p — 2a) pa + paa^), where k is the 
oscillator decay rate scaled by Soj. 

For A ^ 1 the scaled rate of inter-SVS switching 
Wsw/Su! ^ K. Then over dimensionless time ~ there 
is formed a quasi-stationary distribution of the oscillator 
over the states ipn localized around the initially occupied 
SVS. It slowly evolves over time ~ W^sw^- 

From Eq. ([5|, if damping is small, so that n is small 
compared to the transition frequencies |(?n— ffnii |/A, time 
evolution of the populations p„ of quasienergy states is 
described by the balance equation ([T]), where 

W,^:l = 2K{n -h l)|a„„p, W^t""^ = ^W^t- (6) 

n + 1 

Here, a„„ = {i>m\a\ipn) = {2\y^/'^{Q + iP)mn- 

For A ^ 1 one can find a„i„ using the WKB approx- 
imation for ipniQ)- A significant simplification comes 
from the fact that classical trajectories Q{T;g) of the 
system with Hamiltonian g{Q,P) are described by the 
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Jacobi elliptic functions; Q{T;g) is double-periodic on 
the complex-r plane, with real period Tp^\g) and com- 
plex period T^\g) OHl]. For |m — n| <C the matrix 
element is given by the Fourier m ^ n component 
of the function a(T; 5) = [2\)-^/^[Q{t; g) + iP{T; g)]^. 
This gives 



(c) n{n + Vfk'^v^ exp[fci/„ Im (2t» — Tp^-*)] 



f3\ 



I sinh[ifc^'„ /2] \ 



(7) 



Here, t* = T*(g„) and Tp ' = Tp ' (gn) [Im r*,Im ' > 



.(2) 



0]; is the pole of Q(t; f/„) closest to the real axis; 

Vn = 27r/Tp^'' (g„) is the dimensionless frequency of vibra- 
tions in the rotating frame with quasienergy (/„. To the 

leading order in A, we have W^'^+fc = ^n-kn- Equa- 
tion ([7| has to be modified for states very close to the 
extrema of g(Q, 0). 

From Eq. (ITj) , the ratio of the transition rates is 
of the form of a power law, 



_ 2u„ Im (2t,-t 



(2)^ 



(8) 



for I A: I <C A~^. One can show that ^„ < 1. Therefore 
the rates of transitions from a state n toward the SVS 
^^nn+k ■with k <Qi) are larger than away from this state 
(A: > 0). This is to be expected, as approaching the SVS 
corresponds to relaxation. However, even for T = Q there 
are also transitions away from the SVS. 

A power-law transition rate ratio corresponds to de- 
tailed balance. For T = 0, for the quasi-stationary dis- 



tribution we have p^+k/ P' 



(0) 



n-^k n ■ 



III. THE EIKONAL APPROXIMATION 

For arbitrary T we seek the quasi-stationary solution 
of the balance equation ([T]) in the eikonal form 



/9„ = exp(-i?„/A), Rn^R{gn 



(9) 



Away from the critical temperature region R{g) is 
smooth and Rign+k) - R{gn) ~ R'{gn){gn+k - 9n)- 
Then Eq. ([T]) is reduced to a polynomial equation for 
exp[— t'„|_R'(5„)|] [6, 7;. This corresponds to the real- 
time instanton approach in the problem of the distribu- 
tion of reaction systems [16j . The boundary condition is 
R{9n=i) — for A — >■ 0; it corresponds to the occupation 
~ 1 of the closest to the SVS state n = 1. Interestingly, 
\R' {gn=i)\ > 0, in contrast to the conventional instanton 
theory [T7] . 

The dimensionless switching rate M^w is determined 
by the population of states with quasienergies gn close to 
the saddle-point quasienergy gg in Fig. [l] [TH| , 

Wsw ~ Kexp(-i?A/A), 



For T = 
2Im [rj,^\g) 



Ra = Rigs)- (10) 

from Eq. ^ \R'{g)\ = \R"-°H9)\ = 
— 2T*(g)]. This solution is perturbed for 



T > 0. For a given state n, the perturbation is character- 
ized by the ratio >Cn{T) of the rate of absorption-induced 
transitions to this state to the rate of transitions from it. 
For n using the unperturbed distribution we have 



(o) ' 
nm 



z«_ = P^(^t^)pW/^(o). (u) 



From Eq. ([7|), terms Wnm in the denominator in Eq. 
exponentially decay with increasing |n — m|. In contrast, 
Wmn exp [2{n — 'rn)i/nf{gn)] becomes large for n — m ^ 
1 if /(<?„) > 0, 

/(5) = Im [rf (.9)-3n(g)]. 

This leads to the divergence of >f„ for A — > and the 
breakdown of the perturbation theory in n. 

We will analyze the distribution about the small- 
amplitude SVS, see Fig. [l] which is of primary interest 
for the experiment, cf. [101 HI]. In this case f{g) mono- 
tonically decreases with increasing g and goes through 
zero for g = ge\ in our case g^. = 1/4 independent of j3. 
We will consider the range /3 < 2/27, where gs < ge- 

Finding requires evaluating the rates Wn'^m, and 
thus the matrix elements amm for n ~ m > A^^. The 
overlap integral of the wave functions tpn and ipm for 
large |n — to| is exponentially small, which suggests using 
the WKB approximation [TS]. The next step is to use 
the conformal mapping of the complex-Q plane on the 
complex-r plane performed by classical orbits Q{T]g). 
The matrix elements are determined by the behavior of 
these orbits for \Q\ — >■ 00, which allows calculating both 
the exponent and the prefactor in a^n, and thus gives 
Wmn, see Appendix. 

The terms Wmn display a maximum as a function of m 
for gm closest to g^- Calculating J2m ''^mn by the steepest 
descent method, we obtain 



exp 



2A-1 / M)dg' 



(12) 



Cn oc {ge — 9n)'^ is independent of A and T, see Appendix. 

IV. CORRECTION TO THE T = 
DISTRIBUTION AND THE BREAKDOWN OF 
THE INSTANTON APPROXIMATION 

The correction to the distribution can be sought in the 
form i?„ = + A„. For |A„| < A from Eq. Q 

(13) 



-Ax„a„, 



>0, 



where a„ ^ 1. From Eq. (12), |A„|/A increases with n 
exponentially, A„/A„_i « exp[2i^„/(5„)]. The instan- 
ton approximation of smooth R'{gn) ~ ^'^"''(ffn) breaks 
down for |A„|/A ^ 1. Thus the condition >f„ ^ 1, or 
more precisely, jxr^ — 1| is minimal defines the character- 
istic breakdown quasienergy = g^^iT) for given T and 
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the characteristic breakdown temperature = T^{gn) 
for given gn- The thermal perturbation is small for given 
g provided T ^T^{g). 

As T increases it first reaches T^{g) ior g ~ gs. There- 
fore for T < Tc = T^{gs) oc fi? the distribution p„ and 
the switching rate W^sw are described by the T — Q- 
expressions. For T > the T = 0-approximation ap- 
plies only for g — g>({T) A. 

For g between g^tiT) and 175 the distribution is strongly 
changed by the absorption-induced transitions. From 
Eqs. ([1]) and for gn,gm below g^ the corrections to 
R'^^\g) are no longer small, A„ — « — 2At'„/(g„)(n— 
m) for |n — m| ^ 1/A. The transition from the exponen- 
tial dependence of A„ on n for g„ > g^ to this smooth 
dependence occurs in a narrow range around the value 
of n where >£■„ = 1, see Appendix. It corresponds to a 
kink of R'{g) centered at g — g^{T), where R'{g) changes 
from « R'^°\g) to sa i?^(c,) = -2Im T^{g). Such kink 
is indeed seen in the numerical data in Fig. [2] (a). For 
n ^ A'^, absorption- induced transitions to states n with 
g>c — 5n 3> A come primarily from states with g ~ ge, as 
seen from Eq. ( 12 1. 
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FIG. 2. Left panel: the activation exponent Ra calcu- 
lated numerically from Eq. ([T]) (left curve) and the derivative 
dRA/XdT-^ « -dlnWst/dT'^ (right curve). The results re- 
fer to ^ = 0.0448 and A = 0.0053, where the number of local- 
ized states is M = 20. Right panel: the steepness R' — dR/dg 
calculated numerically as —u^^\n{pn/pn-i). The results re- 
fer to the temperatures marked on the left panel by the cor- 
responding symbols. The solid and dashed lines show R' for 
T = (i?'<°') and T -> (R't) in the limit A -s- 0, respectively 

With increasing temperature, the kink of R'{g) at g « 
(7j< moves from gs to ge- From Eq. (12), it approaches ge 
and disappears for n ~ A^/^ or T ^ Tjnst = ^wo/fesl In A|, 
still deep in the quantum domain. For higher tempera- 
tures but still n 1, the slope of R'{g) changes in a 
narrow region \g — gd ^ 1 from i?'(°^((7e) to R'j,{ge)- A 
correction to R'rp{g) for g < g^ can be found in a way 



similar to Ref. 7; it is 



.1/3 



see Appendix. 



From the above arguments, for Tc < T < Tjnst the 
effective switching activation energy Ra ~ A| In Wsw| is 

/g>c ras 
dgR"^''\g)+ / dgR'rrig), (14) 
-svs •'Sx 

(^gyg is the quasienergy of the occupied SVS). It depeds 
on T through g^. From Eq. dlnW^^/dT-^ « 

— (Swo/Zcb) is independent of the modulation. This un- 
expected behavior is confirmed by numerical simulations, 
see Fig. [2] 



V. CONCLUSIONS 

In conclusion, we have studied the probability distri- 
bution and switching of a resonantly modulated non- 
linear quantum oscillator. We find a region of quan- 
tum crossover, which lies between the critical temper- 
ature Tc cx and Tmst oc fi,/|ln/i|. In this region 
the slope of the logarithm of the oscillator distribution 
over quasienergy states, which is the analog of the ef- 
fective reciprocal temperature for this distribution, dis- 
plays a kink. The slope changes from the T = value, 
where the dissipation-induced interstate transitions are 
balanced within a few nearest states, to the value where 
long-range transitions are important. The kink is not de- 
scribed by the conventional instanton theory, which as- 
sumes that the slope of In pn is smooth. As a consequence 
of the kink, in the deeply quantum regime the oscillator 
switching rate is of the activation form with the activa- 
tion energy fejQ independent of the modulation strength. 
The results bear on the current experiments on nonlinear 
cavity modes and modulated Josephson junctions. 
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EMT/QIS 082985, and the ARO, grant W911NF-12-1- 
0235. 



Appendix A: Matrix elements between remote 
quasienergy states 

In this Section we calculate in the WKB approximation 
the exponent and the prefactor of the matrix element 
of the lowering operator between the eigenstates ipm and 
of the quasienergy Hamiltonian g given by Eq. (4) of 
the main text. The approach we propose is similar to, 
but not identical to, that developed [12] in the problem 
of the effect of counter-rotating terms on the dynamics 
of a parametrically modulated quantum oscillator. 

The Hamiltonian function g{Q, P) has the shape of a 
tilted Mexican hat, see Fig. 1 of the main text. For the 
values of g between the local maximum and the saddle 
point of this surface there are two classical phase-space 
orbits P{Q,g), which lie on the inner dome and on the 
external part of g{Q, P). Our analysis refers to the small- 
amplitude stable vibrational state SVS, and we count the 
quantum states off from the state n = 1, which is closest 
in quasienergy to the local maximum of g{Q,P). In the 
spirit of Ref. [TS] 



= Re 



dQa+{Q), 



n > m. 



a+{Q) = Vm(Q)(Q + A9q)V+(Q)/V2A. (A1) 

Function 'ip^{Q) is an eigenfunction of operator g such 
that Re['(/'^((3)] = ipn{Q) and that slightly above the 
real Q-axis in the WKB approximation 



-2i/„ 



1/2 



exp[iX-^Sn{Q) + in/4] (A2) 
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Here, dpgn is dpg calculated for the momentum P = 
PiQjSn), Sn{Q) is the mechanical action counted off 
from the right turning point QR^gn) of the classical orbit 
that lies on the inner dome of g{Q, P) with g{Q, P) = gn- 
For the inner-dome orbits the momentum is 



PiQ,g) = ^jl-Q^-2^g + p^/H 



(A3) 



and5„(g) = /^^(^^)P(Q',g„)dQ'. 

The WKB approximation does not apply to function 
a+((5) close to the zeros and the branching points of 
P{Qi9n) and P{Q,gm)- We go around these points by 
lifting the integration contour above the real Q-axis. 

On the real-Q axis the WKB wave function with 
quasienergy gm is ^pmiQ) c>c cos[iX~^ SmiQ) + iTr/i] for 
Q in the interval between the left and right turning 
points, QL{gm) and QR{gm)- Above this interval on 
the complex-Q plane one of the terms in the cosine be- 
comes exponentially small and should be disregarded in 
the WKB approximation [T5], so that 



i'miQ) 



27rdpgr. 



1/2 



exp [- 



r/4]. 



(A4) 

We will illustrate the calculation by considering the 
case gm < V^- For such g, the orbits on the external 
part of g{Q, P) have the shape of a horseshoe. Function 
P(Qigm) as given by Eq. (A3) has a branching point 
Qb = —g/^fPi which corresponds to the two utmost 
negative-Q points of the horseshoe. It also has three ze- 
ros: Qnigm), QLigm), and Qoxt(5m) on the real-Q axis, 
with Qb < QcKt < Ql <Qr- 



To calculate the integral in Eq. ( Al ) we change in the 
range Q > Qcxt to integration over a contour C shown in 



Fig. [3] (a). On this contour, from Eqs. (All, (A2), and 



(A4) a+(Q) «aWKB(g) ^^^^ 



0^™((5) = TT ^ {vnVm/^XdpgndpgmY^^ 

x[Q + iP{Q, 5„)] exp {i [Sn{Q) - S„,{Q)] /A} . (A5) 

It is straightforward to show using the full expression for 
ipmiQ) that the real part of the exponent in the expres- 
sion for a+((5) monotonically increases on the interval 
(— oo,(5ext]- Therefore, to logarithmic accuracy the up- 
per bound of the contr ibut ion of the integral from —00 to 



Vcxt of a+{Q) in Eq. \kl\ is 
Below we show that 



-(Qoxt) 



,WKB 



(Qoxt)- 



a 



^^^(Qext). Then 



is exponentially larger than 



a-mn ~ Re y dQa 



WKB 



(A6) 



To evaluate the integral ( A6 ) , we analytically continue 



^WKB ^jjg Q-plane with a branch cut on the semi- 
infinite interval {QB{gm),oo). We then can change from 
integration along C to integration along the circle Care 
and contour C shown Fig. [3] (a). 



Classical trajectories Q{t; g) for the Hamiltonian 
g{Q, P) are expressed in terms of the Jacobi elliptic func- 
tions [B]. For each g, this expression provides conformal 
mapping of the Q-plane (with a branch cut) onto a g- 
dependent region on the plane of complex time r. We 
define T{Q,g) as the duration of classical motion from 
the turning point Qpig) to Q. Then the region of the 
T-plane that corresponds to the Q-plane (with a branch 
cut) is the interior of a rectangle shown in Fig. [3] (b). 
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FIG. 3. (a) The contour of integration C for calculating 
the matrix element ( A6 1 in the WKB approximation and 
the auxiliary integration contours C' and Care for > g,n\ 
Qb = Qsigm) and Qcxt = Qext(flm) are the branching point 



and turning points of P{Q, gm), see Eq. (A3 1, (b) Mapping of 
the Q-plane (with a branch cut from Qext to 00, the black thin 
line) on the interior of a rectangle on the r-plane for Qm < 
by function Q(T\g) that describes the classical Hamiltonian 

trajectory with given g, g-n < g < gm', Tp^', Tp^\ and r, are 
the real and imaginary periods and the pole of Q(t; g), respec- 
tively. The solid (C), dashed (Care), and dash-dotted (C') lines 
are the maps of the corresponding contours in (a). The arc 
in the lower left corner is the map of the real axis of Q from 
—00 to Qsig)', tb and Text are the times for reaching QB{gm) 
and Qcxt(ffm)- (c) Integration contours for gm > \f^- (d) 
Conformal mapping for g > y'/J. The curved red line dividing 
the rectangle into two parts is the map of the real axis of Q 
from —00 to QB{g) 



Since T{Q,g) — dS/dg, the exponent in Eq. (A5) is 



dgT{Q,g) . 



(A7) 



As seen in Fig. [s] (b), for any g between gn and gm, 
for any Qarc on contour Care and any Q' on contour C', 
|Im r(Qare,5)l < |Im r(Q',g)|. Therefore, aWKB(Q) 
exponentially smaller on contour C than on contour Care, 
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and the integral along C can be disregarded. Mo reover, 



,WKB 



(QoxtCffm)) < Omn, as assumed in Eq. (A6) 



The integral along Care can be evaluated using the 
asymptotic expressions P « —iQ, gp ~ i2f3^^^Q^^'^ , 



X 



[SniQ) - S^iQ)] 



A^l/4Ql/2 



With Eq. (A8| the integral (A6) is reduced to a simple 



residue, which gives 



/3A3 



1/2 



(5m - 5n) exp 



dgn [g) 



For g„ 



(A8) 



> \fji, the small-momentum branch P{Q,gm) 
(A3) has only two turning points. In the WKB approx- 
imation, a+((5) R:! a^^^{Q) on contour C whose left 
endpoint is the branching point QBigm), see Fig. |3](c). 
The conformal mapping Q{T;g) has a different topology, 
which is shown in Fig.[3](d). Nonetheless, using the same 
arguments as before, we arrive at the same expression 
(A8) for the matrix elements 



In Fig.|4]we compare the explicit analytical expression 
for the matrix elements a,„„, which includes both the ex- 
ponent and the prefactor, with the numerical calculations 
based on solving the Schrodinger equation gipn — gnipn- 
The results are in excellent agreement. 




0.22 0.24 

9n 



term Wmn = Wmn^^ Pm^ / Pn ^ in the expression for >f„ for 
n — m 3> 1 is of the form 



f{g)dg 



C^'^\g„.,gn) 



(Bl) 



As a function of to, w„in is maximal for gm closest to the 
quasienergy value ge given by the condition f{ge) — 0. 

The leading-order contribution to >Cn comes from the 
terms w,nn with \gm — gd *C 1- The sum of these terms 
can be calculated by changing to integration over m 
and then using the steepest descent method. This gives 
Eq. (12) of the main text for >r„, with 

C„ = C^-\g,,gr.)[K\/f{g,)]'/'/ [xwi'^) , 

W^/^^^E ^nrl (B2) 

Here, /' — df /dg, whereas Wn^ is the overall rate of 
transitions from state n for fi <C 1. It can be found from 
Eq. (7) of the main text; using this equation we obtain 

C„^CW(ge,g„)k//'(5e)]^/^ ^ 

We now consider the distribution for T > 0, 



exp(-i?„/A), 



(B3) 



For small n we need to consider the absorption-induced 
transitions to a given state n only from states to with 
m < n. Then the balance equation in the quasi- 
stationary regime reads 



= E ^ ( 



c) \ g(A„-A„.)/A 



^/w, 



1 = -Pn 
(A„-A„0/A^ 

0. 



(B4) 



The parameter >f„ is given by p„ for A 



1. Temperature range n <C 



FIG. 4. Comparison of Eq. (A8l for amn calculated as a 



continuous function of Qn (solid lines) with numerical calcula- 
tions (symbols). The scaled intensity of the modulating field 
is P = 0.0448. The other parameter values are A = 0.0053 and 
m = 4 (violet triangles), A = 0.09 and m = 2 (blue circles), 
and A = 0.15 and m — 1 (green squares). Parameter m has 
been chosen so that ~ (?e = 1/4. 



Appendix B: Thermally-induced modification of the 
T — distribution 



The explicit expression for the transition matrix el- 
ements (AS) makes it possible to calculate the rate of 



absorption-induced transitions to a state n scaled by the 
rate of leaving this state Kn{T), Eq. (11) of the main 
text. The most probable transitions are those from states 
TO closer to the SVS, m < n. 



From Eq. (AS I, the 



Thermal modification of the distribution becomes sub- 
stantial when n is exponentially small in A, |lnn| > 
R^^\gs)/\. In this subsection we consider the distri- 
bution for fi, <C A'^. 

We start with the range of the level numbers n where 
n^ — n^ 1; the quasienergy level number n^^ is given by 
the condition that \xn^ — 1| be minimal, gn^ = gx{T). 
For such n, |A„|/A ^ 1, and in p„ one can disregard 
An, Am- Then p„ — Xn is given by Eq. ( 12) of the 
main text. In the left-hand side of Eq. (B4) it suffices 
to keep the terms linear in A„, A,„. Taking into account 
that Xn/xn-i ~ exp[2i/„/((7„)], from Eq. (B4) we then 
obtain Eq. (13) of the main text with 



1 



•i)/(9,0 



Using the explicit form of Wnm and taking into account 

f2) f2) 

that Im {jp ' — 2t*) > |Im (t^ — 4r*)|, one can show 
that a„ > 0; clearly, a„ is independent of A and n. 
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We then consider the states with n > n^, and first 
assume that n — n^ <^ 1/A. In this range it is convenient 
to spht the sum over m in p„ in Eq. ( B4 ) into a sum from 
n — 1 to rij; — no and a sum over m < — hq, where hq is 
chosen so that I ^ no ^ 1/A. Since A™ is exponentially 
small for rt^ — m 3> 1, it can be dropped in the terms 
with m < rijj — tiq. Taking into account the explicit form 
of Wmn, Eq- (Bll, and that, as a consequence of this 
equation, >«:„ depends on n exponentially, we write p„ as 



Pn ^ 1 



,(A„-A„)/A 



g(A„/A)+2(n-n^)i.„/(g„)^ 



From the relation >c„ 



1, it follows that the solution of 



Eq. (B4) is 



A„ - A„i « -2Xvn{n - m)J{g.n),n- m < A \ 
A„ « -2A(n - (B5) 



provided n — n^ > m — n^ 3> 1. Corrections to A„ come 
from the terms A„' with n' — ~ 1, |A„'| ~ A. They 
also come from the prefactor in Wmm the dependence of 
Vn and Qn on n, and the dependence on n of the left-hand 
side of Eq. ( B4 ) calculated for A„ — A^ of the form ( B5 1 . 
All these correction are ~ A. For A„ of the form (B5) 



P 



(1) 



< 



(2) 



In contrast to the exponentially steep n-dependence of 
A„ = — Aa„>f„ for n < n^, in Eq. ( |B5| ) A„ smoothly 
depends on n. The region of n where the two expressions 
join one another is centered at n^. The width of this 



region is independent of A or n, as seen from Eq. ( B4 1 
In this region A„/A ^ 1. 



1/A. Here, too, pn ^ < 



One can easily see that Eq. (|B5|) applies also for n 



n{n ■ 



l/X^ < p'n 



(2). 



in the expression for A„ one should replace 



2{n ~ n^)v,J{gn) 2A" 



dgfig). 



The inequality pn"^ <^ pn^ indicates that in the tem- 
perature range n <C A'^ absorption-induced transitions to 



states n > come primarily from remote states with 
quasienergy w g^- 

2. Temperature range n ^ A^ 

As the temperature increases n becomes larger than 
A'^ and the term pn'' becomes more important. Equa- 
tion ( |B5| still gives the leading-order term in (A„ — 
Am)/A in the range n — n^,m — n^^ ^ ttT^I"^. How- 
ever, when calculating p^n^ one should add a correction 
— (n — m)e„ to this term, |e„| ^ 1. From Eq. (Bl), the 



sum over m in p^n^ is then a second derivative with re- 
spect to e„ of a geometric series, since in C^"^^ we have 



(.9n 



9r, 



\\vnin — m)]2. The summation gives 



Pi'^ « (^/epC W , C« = 8^^/A/3I^(<=) . (B6) 
The coefficient Cpn is independent of A and is ^ 1, gen- 
erally. Since we need pn < I, from Eq. ( |B6[ ) we have 
e„ ~ Formally, the sum over rt — to in pn^ was 

extended to infinity, which requires that n - 



and then from Eq. ( B5 ) e„ ^ A 



> 1/en, 



(2) 

The correction to A„/A in the expression for p„ is 
given by J2k<n^n- It becomes ^ — n~^/^/A 3> 1 for 

n — n^ ^ A^^. Therefore pi^"* is small for such n. On the 
whole, as the ratio n^/^/A changes from small to large, 
so does also the ratio pn \ starting first from large 
n and then going to smaller and smaller n. 

From Eq. (Bl ) and from Eq. (12) of the main text, for 
fi approaching A'^/^ ^ i the quasienergy g^ approaches 
ge- For n ^ \^^^ (T > Tingt) the rates of transitions 
into all states with g„ < g^ are determined by thermal 
processes. In this range p„ = pi^"* and the m odifi cation 
of Rn compared to is determined by Eq. (B5| 

Keeping in mind that i?„ — Rn-i 
from n « 



wc obtain from Eqs. (B3) and (B5) 
R'{g) « R'^g) 



-XR'ignJpn away 
(B7) 



-2Imn(g), g < g. 



Thus in the whole transition region where 175 < < 
5e (and respectively, Tc < T < Tinst) the distribution 
depends on temperature primarily through the position 
of the kink of R'{g), where R' sharply changes from 
to R'. 
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